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Abstract 

The phenomenon of stochastic resonance (SR) is known to occur mostly in bistable systems. 
However, the question of occurrence of SR in periodic potential systems is not conclusively resolved. 
Our present numerical work shows that the periodic potential system indeed exhibits SR in the 
high frequency regime, where the linear response theory yields maximum frequency dependent 
mobility as a function of noise strength. The existence of two (and only two) distinct dynamical 
states of trajectories in this moderately feebly-damped periodically driven noisy periodic potential 
system plays an important role in the occurrence of SR. 
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I. INTRODUCTION 



Stochastic resonance (SR) was discovered theoretically about three decades ago[l|. Since 



then SR has 
instance refer 



3een investigated with vigour and many notable reviews have appeared, for 
2j, y|. Physical systems are always subject to internal or external (thermal or 
otherwise) fluctuations (noise). The optimal periodic response of a system to an external 
periodic drive as a function of noise strength is termed as stochastic resonance. It has been 
experimentally found to occur, to just mention a few, in electronic circuits {4 -6], two-mode 
ring lasers^, nanomechanical systems)^, neuronal systems {9- 12], etc. Its main attraction 
lies in its practical utility of selecting and enhancing a signal of a particular frequency out 
of a host of signals by tuning the noise strength. Presumably, biological systems use SR to 
their advantage 13j . It has the potential to be utilized to control kinetically the pathways 



of a biochemical reaction 

SR has been predicted and shown to occur mostly in bistable systems^, flo| . 



However, 



there has been some notable investigations of SR in monostable and also periodic potential 
systems [If], \v\ . In the monostable systems SR is shown to occur in the high frequency 
regime close to the natural frequency of oscillation at the bottom of the potential. However, 
the occurrence of SR in periodic potentials have not been conclusive [18]. 



Dykman and coworkers 



171 ] use an interesting model in which a single-well Duffing oscilla- 



tor with additive noise is driven at a frequency close to the natural frequency of the oscillator. 
The model is used under linear response theory formalism to study fluctuation phenomena 
associated with two coexisting periodic attractors. A weak Gaussian noise causes transition 
between these two attractors. The populations of these two attractors become equal at a 
particular noise strength where the response becomes maximum. This is considered as a 
genuine signature of a (nonconventional) stochastic resonance. The theoretical result was 
supported by an analog electronic circuit experiment. A similar result was obtained in an 
underdamped superconducting quantum interference device fiol]. 

The same resonance behaviour, as in Ref.[l7|], in the frequency dependent mobility was 



obtained in a periodic potential using linear response theory by Kim and Sung[18| in the high 
frequency range of the external periodic drive. However, these authors ascribe this resonant 
behaviour as simply a noise assisted standard dynamical resonance as the transitions involve 
only intra-well motion. Also, in the inter-well hopping (low) frequency regime the frequency 
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dependent mobility shows monotonic behaviour as a function of noise strength thereby 
discounting the possibility of occurrence of SR in periodic structures. However, the authors 
show that SR can occur if the driven system has a tilted periodic potential so that the 
passages are allowed only in one direction. 

Moreover, recently, it was observed that in a bistable potential, V(x) = V e~ ax2 + b\x\ g /q, 
the confinement parameter q plays an important role in deciding whether the system will 
show SR or not 20]. For q > 2 the system shows SR whereas for q < 2 it does not. In 
addition, we find the input energy expended per period of the external field on the system 



by the field acts as a good quantifier of SR 



211 ] . This input energy is ultimately dissipated 



into the thermal bath. This is naturally a measure of the hysteresis loop area in positions- 
force^) space. Although the input energy and hysteresis loop area are exactly the same 
in magnitude the latter is an average quantity, whereas input energy has a well defined 
distribution. The input energy distribution provides useful information about stochastic 
resonance behaviour. In particular, the distribution shows a characteristic largest shoulder 
(bimodality) at stochastic resonance 22h24| . On the other hand, hysteresis loops carry 
important information about phase relationship between x and F which have also been of 



interest to SR 



21 



25 



26]. 



In the present work, we explore the possibility of occurrence of SR in a periodic sinusoidal 
potential using input energy and hysteresis loop quantifiers. Moreover, one can 

take various values of the wavevector k of the potential analogously varying the effect of 
confinement parameter q of the bistable potential discussed above. We, however, present 
results for k = 1 only. 

We find that the periodic sinusoidal potential does not show SR when driven by a low 
frequency field corresponding to Kramers rates across the maxima of the potential or when 
driven at a still lower frequency. The same conclusion have been arrived at in Ref. 27] 
while studying the diffusion coefficient in a periodic system. However, it should be noted 
that in Ref. 27j the probability -P(t') of a particle, after going from one well to an adjacent 
one returns back to the same initial well in the subsequent time r', shows periodic peaks. 
The strength of the first peak of P(t') shows SR-like behaviour. However, -P(r') has, by 
construction, the bearings of a bistable system and not that of a periodic potential system. 

We further find that in the high frequency range the input ener gy behaves exactly simi- 
larly as the response function in the works of Dykman and coworkers 17j and as the frequency 
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dependent mobility does in the work of Kim and Sung 18j. In addition, our work shows that 
the input energy peaks as a function of noise strength. This is an indication of SR arising 
due to a competition between two dynamical states of particle trajectories (to be elaborated 
in section III) as in the case of bistable systems. Though the trajectories are intra-well 
in nature close to SR the transition between these two states are also aided by inter-well 
passages of particles across the potential maxima. 

The two dynamical states of trajectories are distinctly identified by the phase difference 

between the periodic forcing F = F cos(cot) and the trajectory x(t) = xo cos(ut + 0); one 
having a fixed phase lag <p — <pi and the other = 2 - Note that the system, at finite T, 
being stochastic in nature and 02 are average quantities. These individual phase lags 

01 and 02 effectively do not vary with the noise strength. However, the relative cumulative 
length of these two dynamical states in a trajectory change with the noise strength. The 



average phase lag, there 



and observed in Ref. 



25 



ore, changes with noise strength, similar to what was predicted 



26[ in the case of bistable systems. Moreover, the distribution of 
input energy shows very similar behaviour across the input energy peak as in case of SR 
in a bistable system thus affirming the genuineness of SR in the present periodic potential 
system. 

We consider two model systems for our study. In one case the medium is considered to 
having uniform friction, whereas in the other case the friction is considered to be nonuniform. 
The nonuniform friction 7(2) has an exact analogy in the Resistively Coupled Shunted 
Junction (RCSJ) model of Josephson junction. In the RCSJ model the 'cos0' term, arising 
out of the interference between the quasiparticle tunneling and the Cooper pair tunneling 
across the junction, is equivalent to the nonuniform part of j(x) in the present case. The 
systems with nonuniform friction, however, are not very rare and are not without practical 
relevance [28j|. Since the potential is symmetric and periodic the homogeneous system does 
not show any average mobility. The nonuniform system, however, shows average current 
when driven by a sinusoidal forcing in the underdamped system 29]. 

Interestingly, for the uniform system, 0i ~ — 0.0137T and 02 — — 0.57T whereas in the 
nonuniform case, 0i ~ — 0.0257T and 2 ~ — 0.857T. In either case the trajectories have dis- 
tributions of these two phases depending on the noise strength. This is reflected in the form 
of the (x — F) hysteresis loops and hence in the average response amplitude and phase lag. In 
the homogeneous case the behaviour of the average phase lag only approximately conforms 
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io the SR prediction of 2l|, |26|] and in the other case it follows closely the observations of 
25[ in bistable systems. 



II. THE MODEL 

In this work, we consider the underdamped motion of a particle in a periodic potential 
V(x) = —Vq sin(fcx) which is symmetric in space (about kx = (2n+l)7r/2, n = 0, ±1, ±2, ...). 
The system is driven periodically by an external forcing F(t)=F cos(ut). We study two cases 
of the system - when the friction coefficient 7(2) is uniform (=7o) (system is homogeneous) 
and when the friction coefficient is space dependent 7(2) = 7 (1 — Asin(/cc + 9)) (system 
is inhomogeneous) . In the latter case, the friction is periodic with the same periodicity as 
the potential but has a phase difference 9 with it (9 ^ 0, 7r). A (0 < A < 1) determines the 
degree of inhomogeneity of the system (A = corresponds to the homogeneous system). 

A particle of mass m moving in a periodic potential V(x) = — Vq sm(kx) in a medium 
with friction coefficient 7(2) and subjected to an external periodic forcing F(t) is considered 
to be described by the Langevin equation, 



d 2 x dx dV(x) 



m—— = —j(x) 



dt 2 x dt dx 

d 2 x dx dV (x) 

m ~iw = _ 7o- 



+ F(t) + v ^jT^(t). (2.1) 



dt 2 -dt ox +nt) + yfirm (2.2) 

Eqn. (2.1) is for the inhomogeneous system and Eqn. (2.2) for the homogeneous system. The 
temperature T is in units of the Boltzmann constant ks- The inherent random fluctuations 
in the system are represented by £(t) which satisfy the statistics: < £(t) >= 0, and < 
f (t)f (i ) >= 25 (i -t). The equations are written in dimensionless units by setting m = 1, 
Vo = 1, k — 1. The Langevin equation, with reduced variables denoted again now by the 
same symbols, corresponding to Eqns. (2.1) and (2.2) are written as 

d£ = ~ 7(X) + C ° SX + F(t) + VtWTZV)- (2-3) 
72 



cf/ oc doc 

= ~ l0 ltt + C ° SX + F(t) + V^^W- (2-4) 
The potential barrier between any two consecutive wells of V(x) disappears at the critical 
field value Fq = F c = 1. The noise variable, in the same symbol £, satisfies exactly similar 
statistics as earlier. 
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III. NUMERICAL RESULTS 



The trajectories x(t) are obtained numerically 32j by solving the Langevin equations (2.3) 
and (2.4) corresponding to the inhomogeneous and the homogeneous system, respectively, 
with the amplitude of the drive field, F Q = 0.2 and u> = 2tt/t, with r = 8. The value of 
oj is close to the natural frequency of the potential. At high temperatures T the average 
behaviour of these trajectories is the same independent of the initial conditions. However, at 
low temperatures, and especially in the limit of deterministic motion, and at such a low field 
amplitude as F = 0.2, the trajectories are intra-well in nature. Yet, their behaviour is very 
sensitive to the initial conditions, x(0) = x(t = 0) and v(0) = v(t = 0) 30] . In our numerical 
calculations we take v (0) = and x(0) at N equispaced intervals, x^, i = 1, 2, N, between 
the two consecutive peaks, e.g. [— n/2 < Xi < Stt/2], of the periodic potential V(x). In most 
of the cases we take N = 100, but at some values of temperature we take N = 300. 




-°- 5 400 



420 440 460 480 500 520 540 
t 



FIG. 1: Plot of x(t) (solid line) and F(t) (dot-dashed line) for T = 0.015. The particle exhibits 
both kinds of trajectories with a transition from out-of-phase to in-phase, in this case, at around 
t = 460; r = 8, F = 0.2, 7o = 0.12. 



Depending on the range of x(0) we get two distinctly different kinds of trajectories; 
one that lags behind the applied field by a small phase difference <j>\ and the other by 
a large phase difference 02- Only for the sake of convenience, we call the former kind 
of trajectories as being in-phase and the other as being out-of-phase. The out-of-phase 
trajectories always have much higher amplitude than the in-phase trajectories, Fig.l. At 
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the lowest temperature considered here there is no transition between these dynamical states 
of trajectories. However, as the temperature (or the noise strength) is increased transitions 
do take place, Fig.l. Yet, the trajectories are found basically only in these two states even 
at temperatures where inter-well transitions lead the trajectories out of the initial well. 

Following stochastic enegetics formulation of Sekimoto j3l| the input energy, or work done 
by the field on the system, W, in a period r, is calculated as 

w(toM+T) =r°jiM}A dt , (3 ,) 



at 

where, the potential U(x(t),t) = V(x) — xF(t), and V(x) = — sin(x), F(t) = F cos(cut). 
The average input energy per period over an entire trajectory W, is 

. n=N\ 

W= Af E W(nr,(n + l)r). (3.2) 

1 n=0 

Typically, the number of periods, Aq, taken in a trajectory ranges between 10 5 to 10 7 , as 
required. Finally, the average input energy per period < W > is calculated by averaging 
W over all the trajectories. From Eqn. (3.1) the distribution P(W) is also calculated. 
As we shall see below P(W) provides an important complementary criterion for stochastic 
resonance. 

Eqn. (3.1) can be further written as 

W(t , t + r)= f° +T dU ^^ dt = - [ F{t ° +T) xd F, (3.3) 

Jto 0t JF{t Q ) 

which is the hysteresis loop area over a period. Since x(t) is a stochastic variable it is 
unresonable to expect a sensible hysteresis loop over a period of the field. However, when 
averaged over the entire duration of a trajectory a well defined hysteresis loop x(F(ti)) and 
its area A is obtained: 

j n=N 1 

x(F(t i )) = w J2 x ( F ( nT + t ^ ( 3 - 4 ) 

1 n=0 

for all [0 < ti < t] and 

A=\W\. (3.5) 

The calculation of the hysteresis loops x(F(ti)), Eqn. (3.3), is correct at the lowest tem- 
peratures where the trajectories maintain the same phase <fi throughout and also at higher 
temperatures where trajectories change phase between <fti and 02 during their journey. How- 
ever, the hysteresis loops will be different for different trajectories depending on the cumu- 
lative duration of the dynamical states 0i 2 of the segments of the individual trajectories. 



Therefore, an ensemble average < x(F(ti)) > is taken over all the trajectories considered. 
The hysteresis loop < x(F(ti)) >, with area < A >, reflects the average response x(t) 
of the system to the applied field F(t). In the linear response regime one would expect 
x(t) = xq cos(cut + (p) for F{t) = F cos(ut), u = 2ii/t, where is the average phase lag to 
be observed experimentally |25l] . Notice that will, in general, be different from the fixed 
phases 0i and 02 characterising the individual dynamical phases of the trajectories. 

In the following, we present the results of our numerical calculations and analyse how the 
particle trajectories in the sinusoidal potential and hence the input energy, and the hysteresis 
loops are affected by noise strength. We examine the occurrence of stochastic resonance in 
the periodic potential, referring to the standard SR criteria in bistable systems. 

A. Homogeneous Systems 

1. The intra /inter- well transitions and input energy 

Fig. 2a shows the plot of input energy W averaged over a trajectory with initial position 
x(0) at the lowest temperature T = 0.001 for the homogeneous system (7(2;) = 70) • W 
are confined to two narrow bands around 0.094 and 1.179 corresponding to x(0) values that 
lie in contiguous regions in the range [— 7r/2, 3n/2] and no points in between. The input 
energies W = 0.094 correspond to in-phase trajectories with phase lag 0i ~ — 0.0137T and 
those with W = 1.179 correspond to out-of-phase trajectories with phase lag 02 — — 0.57T. 
The trajectories continue to be in the same state of phase lag throughout their course of 
journey. At this low temperature the trajectories are quite stable. The inset of Fig.2a shows 
the (x — v) plane Poincare (stroboscopic) plots revealing the two attractors corresponding 
to the two dynamical states of the trajectories. The (stochastic) spread in the plots is due 
to the finite temperature of the system. 

It may be recalled that the results are obtained by taking v(0) = 0. However, the two 
dynamical states of trajectory are genuinely stable as revealed, in Fig.3, by the basins of 
attraction of the two attractors at T = 0.001. In the figure the in-phase and the out-of-phase 
states are indicated respectively by inp(Z) and outp(m) for /, m = 0, ±1, ±2, etc. The indices 
I, m (within the braces) specify the well number where the particle settles down as it begins 
from the initial well (l,m = 0). Even though the attractors correspond to wells l,m ^ 0, 
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FIG. 2: Plot of W with x(0) for T = 0.001 (a), T = 0.003 (b) and T = 0.016 (c), r = 8, F = 0.2, 
7o = 0.12. The insets in Figs, (a) and (c) show the stroboscopic (Poincare plots in the (x — v) 
plane at times r n = nr, n = 0, ±1, ±2, ... . 
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the phase and amplitude of the trajectories remain exactly the same as when in the initial 
well (l,m = 0). Thus, the differently labeled trajectories are physically the same two and 
the only two: the in-phase and the out-of-phase with fixed <f>\ and 02, respectively. 

Fig. 2b shows the same plot as Fig. 2a but at T = 0.003. The figure clearly shows that 
a few points have deserted the upper band of W and all of them but one have already 
joined and the lone one on its way to join the lower band of W. This shows that all of 
these few out-of-phase trajectories have made a transition to the in-phase state (as shown 
in Fig.l) quite early during the 10 5 periods of the trajectory and the lone one somewhere 
in the middle of all these periods. It is to be noted that no transition has taken place 
from the in-phase state to the out-of-phase state. This shows that the out-of-phase state is 
less stable than the in-phase state and are separated by an energy barrier of about 0.003 
(in units of Vo, the amplitude of the sinusoidal potential) from the out-of-phase side. This 
trend continues till T = 0.016 where W consists of only one band, that is the lower band, 
Fig.2c. The stroboscopic plot in the (x — v) space, inset of Fig. 2c, shows the lone attractor at 
T = 0.016. The other attractor has now ceased to exist. From T = 0.003 to T = 0.016 the 
ensemble averaged input energy < W > decreases very rapidly attaining the lowest value at 
T = 0.016. 

For T > 0.016, transitions start taking place from the in-phase to the out-of-phase state 
too. Consequently, < W > begins rising for T > 0.016. Of course, the out-of-phase state 
lives for a very short duration before a larger temperature is reached. This is demonstrated 
in Fig.4, for T = 0.04. Therefore, there is an energy barrier of roughly about 0.016 from the 
in-phase state side to the out-of-phase side. One can thus roughly pictorise the two states 
as shown in Fig.5 with the bottom of the wells at —<\>\ = 0.0137T and —02 == -57r. 

As the temperature is increased further the relative population of out-of-phase state keeps 
increasing with increase of temperature. The transitions in both the directions maintain a 
constant ratio at any given temperature. Moreover, by the temperature T ~ 0.08 the inter- 
well transitions have already set in and have become numerous. These inter-well transitions 
also help intra-well transitions, Fig. 6. At around T = 0.2, the relative population of both 
the states become almost equal and the input energy < W > peaks. T = 0.2, thus, falls in 
the region of kinetic phase transition 33J| between the two dynamical states. 

Beyond T = 0.2, the intra-well transitions become more frequent and so do the inter- 
well transitions. The transition points or the attempted transition points lead to lower 
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x(0) 

FIG. 3: The Basins of attraction of the in-phase inp and out-of-phase outp attractors at temper- 
ature T = 0.001 are shown. The bracketed numbers on inp and outp indicate the well number 
of the periodic potential where the trajectory settles down in the in-phase or out-of-phase states. 
The region marked 'a' correspond to inp(-l), 'b' correspond to outp(-2), 'c' to the continuation of 
inp(-2) indicated on the left side, 'd' to inp(-3), and the region 'e' corresponds to outp(-3). The 
boundaries correspond to the end of the phase indicated above the respective lines. The thick line 
separating outp(2) and inp(l) indicates that some other phases too chip in between. The phase 
inp(l) at the top right corner is just the continuation of the phase inp(l) on the left top corner of 
the figure. For this figure, r = 8, F = 0.2, 70 = 0.12. 

amplitude motion. This results in gradual lowering of < W > with temperature. At much 
higher temperatures T ^> 0.2 the (intra-well as well as inter-well) transitions become so 
numerous that the phases effectively lose their distinct identity. This gets reflected in the 
input energy distribution P{W) as a single peak (Fig. 7, T = 1.0). 

In Fig. 7 the probability distribution of W is drawn for various temperatures. At the lowest 
temperature T = 0.003 we see two distinct peaks of P(W): The low W peak corresponding 
predominantly to the in-phase state of trajectory and the other to the out-of-phase state. As 
the temperature is gradually increased the large W peak shrinks and at T = 0.016 this peak 
disappears completely. However, as the temperature is increased further the out-of-phase 
peak reappears and begin to swell to maximise at T ~ 0.2. At T = 0.2 the broad shoulder 
of P(W) characterises the peak in < W >. At the largest temperature shown T = 1.0 both 
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FIG. 4: Plot of the particle trajectories x(t) for T = 0.04. Transistions are seen between the 
in-phase to the out-of-phase states. Inset is a magnified figure. F(t) (dashed line) is also included 
for comparison; xO = 2.0, r = 8, Fq = 0.2, 70 = 0.12. 




0.01371; -<)> 0.57C 

FIG. 5: A schematic depiction of the two dynamical energy states of the particle trajectories. In 

the figure E3 - El ~ 0.016 and E3 - E2 ~ 0.003. 



the peaks merge into a broad single peak. Though < W > always remains > thus never 
violating the second law of thermodynamics, interestingly, P(W) is not confined to W > 
but a significant part of it lies at W < 0. 

Fig. 8 shows the variation of < W > as a function of temperature T. < W > peaks at 
a temperature around T = 0.2. This is a clear signature of stochastic resonance even if 



the criterion for SR in a bistable system is adhered to as a benchmark 



20 



22 



23|. This is 



supported by the behaviour of P(W) across T = 0.2. P(W) shows a prominent shoulder 22 



12 




x 5 - 







1000 2000 3000 4000 5000 6000 7000 8000 
t 



FIG. 6: Plot of x(t) for T = 0.08. The figure shows inter-well transistions as also numerous 
transistions between the in-phase (lower amplitude) and out-of-phase (higher amplitude) states. 
At around t = 3700 and t = 5400, the inter-well jump also leads to a transistion from the out-of- 
phase to the in-phase state; x(0) = 2.0, r = 8, Fq = 0.2, 70 = 0.12. 
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(a) 

FIG. 7: Plot of P{W) for different values of T (Fig. a for T < 0.016 and Fig. b for T > 0.016) for 
the homogeneous system. For large T, e.g., at T = 1.0, the distribution has a single peak structure; 
r = 8, F = 0.2, 70 = 0.12. 



231 ]. characteristic of SR, around T = 0.2, Fig. 7. In the following we examine the behaviour 
of phase lag of the response to the periodic field. 
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FIG. 8: Plot of (W) as a function of T for the homogeneous system; r = 8, Fq = 0.2, 70 = 0.12. 



2. The Hysteresis loop: Area and the phase 



From Eqns. (3.3), (3.4) and (3.5) the hysteresis loop area < A > being equal in magnitude 
to < W >, < A > does not provide any additional information than < W >. However, the 
hysteresis loop itself does give an important insight. The amplitude F = 0.2 of the external 
applied field F(t) = F cos(ut) being small, one expects the response, namely the position 
x(t), to have a linear variation with F(t): x(t) = Xq cos(ut + 0). Hence, the hysteresis loop 

< x(F(ti)) > will closely be an ellipse. The average phase difference is measured from the 
resulting ellipse, Eqn. (3.4). Unlike <pi ~ — 0.0137T, and 02 — — 0.57T, which remain more or 
less constant with T, the average phase difference <p varies with T, as shown in Fig.9. In 
Fig.9, the average amplitude Xq and the average area < A > are also plotted for comparison. 

At the lowest temperature T = 0.001, <fi is close to 02 — — 0.57T. At this temperature, Xq 
as well as < A > are also at their respective maxima. As T is increased —0 as well as < A > 
decrease sharply, and attain a minimum at around T = 0.016. The minimum is close 
to 0i ~ — 0.0137T. This is because at T = 0.016 all the trajectories are in the dynamical 
state of phase 0i. Thereafter for T > 0.016, —0 increases monotonically. The inflection 
point of — 0(T) occurs at T ~ 0.09, where —0 ~ tt/4 as observed in Ref[25]. However, 

< A > becomes maximum only at a much higher temperature T ~ 0.2, where —0 ~ 0.357T. 
Therefore, < A > here does not exactly satisfy the additional approximate SR criterion 
on suggested in Ref 25] . However, in the inhomogeneous system, the SR criterion on 
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FIG. 9: Plot of (A) (top), xq (middle) and —(f) (bottom) with T for the homogeneous system; 
t = 8, Fq = 0.2, 70 = 0.12. Beyond T = 0.4, the curves have a monotonic behaviour and are not 
shown. The top figure is an exact reproduction of Fig.8. 
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suggested using linear response theory 



261 ] appears to be respected. 



B. Inhomogeneous Systems 

In this case the particle experiences a nonuniform friction as it moves in the medium. 
As stated earlier the potential is considered periodic: V(x) = Vosin(fcr). The friction 
coefficient 7(2) is taken as j(x) = 70(1 — Asin(/cx + 6)) instead of the constant friction 
coefficient 7(2;) = 70 as in the homogeneous case. We take fixed values A = 0.9 and phase 
difference 9 = 0.35 throughout. 0, n) provides the necessary asymmetry in the system 



to yield an average particle current, the ratchet current [3 0|, even when driven by a zero 
time-average external forcing F(t) = F cos(cot). For k = 1, 70 = 0.12, F = 0.7, and 
T = 0.4 the system shows a current maximum at u — 27r/140. Keeping u = 2tt/t,t = 140 
and other parameters fixed the current maximizes at T ~ 0.25. However, we do not see any 
maximum in the periodic response, such as the hysteresis loop area, to the external periodic 
field corresponding to these parameters. The system does not show SR at this low frequency 
of drive. 

Here, the particle position x(t) does roughly follow the periodic field with an irregular but 
large amplitude. There is also a small amplitude high frequency component superimposed 
on the low frequency response to the field. The frequency of the superimposed component 
is close to the natural frequency of the potential. In the following we consider the external 
field F(t) with Fq = 0.2 and r — 8, exactly as in the homogeneous case, but keeping in mind 
that in the inhomogeneous case the system is no longer symmetric. Therefore, even at this 
high frequency the system shows an average ratchet current but the current itself does not 
show any peaking behaviour with T. 



1. The intra/ inter-well transitions and input energy 

As in the homogeneous case, the system shows two distict dynamical states of trajectories. 
One with phase difference <fii ~ — 0.0257T and the other with 2 ~ — 0.857T. We again call the 
trajectories with phase difference <\>\ as being in-phase and the other as out-of-phase with 
the external field. The latter having much higher response amplitude xq than the former. 
The frictional nonuniformity, surprisingly, leads to out-of-phase trajectory amplitude about 
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three times larger than in the homogeneous case and therefore larger average energy < W >. 
As a consequence, the results are qualitatively different from what was observed in the case 
of uniform friction system. 
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FIG. 10: Plot of W with x(0) for T = 0.001 (a), T = 0.012 (b), T = 0.014 (d) and T = 0.016 
(d), r = 8, Fo = 0.2, 70 = 0.12 for the inhomogeneous system. The insert of Fig.a shows the two 
attractors in the (x — v) plane stroboscopic plots at times nr. 

At the lowest temperature T = 0.001 the trajectories are bunched into two contiguous 
groups and continue to be so for all time without ever jumping across the groups, Fig. 10a. 
The inset of Fig. 10a shows, in the (x — v) plane stroboscopic plots, the two attractors 
corresponding to the two dynamical states of particle trajectories. The situation remains 
the same till T = 0.01. However, the average energy < W > decreases with T due to 
a slight decrease in the out-of-phase trajectory amplitude. Contrast this with the sharp 
decline in energy in the homogeneous case due to transition from out-of-phase states to the 
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FIG. 11: Plot of < W > (top), the amplitude xo (middle) and the phase difference — (bottom) 
with T for the inhomogeneous system; r = 8, Fq = 0.2, 70 = 0.12. 

in-phase states. At T around 0.011 transitions begin to take place from the in-phase states 
to the out-of-phase states, Fig. 10b; the out-of-phase state is more stable. The < W > rises 
sharply thereafter. However, by T = 0.014 transitions are also observed from out-of-phase 
states to the in-phase states, Fig. 10c. Therefore, the potential barrier from out-of-phase 
side to the in-phase side is about 0.011, whereas in the reverse direction the barrier height 
is about 0.014. Both the states are almost equally stable. By the temperature T = 0.016 
even inter- well transitions are observed, Fig.lOd. 

As the temperature is increased to 0.018 the average energy < W > attains a maximum. 
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FIG. 12: Plot of P(W P ) for different values of T (Fig. a for low T and Fig.b for higher T) for the 
inhomogeneous system; r = 8, Fq = 0.2, 70 = 0.12. 

At this point the two states become almost equally populated, a point in the region of 
kinetic phase transition [33|. The variation of < W > is plotted in Fig. 11 as a function 
of temperature. This is a signature of stochastic resonance occurring at T = 0.018. The 
distribution P(W) again shows the largest asymmetry at the same temperature T = 0.018. 
Here P(W) is distinctly bimodal in nature, Fig. 12. For this nonuniform friction case, P(W) 
behaves in a manner similar to the case of homogeneous system (Fig. 7). However, here the 
out of phase state has the lowest energy unlike in the homogeneous case. At T = 0.018 we 
have the largest contributing second peak of P(W) to < W > at any nonzero temperature. 
These intra-well transitions responsible for SR are effectively supported by numerous inter- 
well transitions in the periodic potential system. 



2. The Hysteresis loop: Area and the phase 

Since the system is asymmetric and the amplitude Xq are large, at low temperatures the 
hysteresis loops are not perfectly elliptical. Yet, it is possible to roughly calculate the average 
amplitude xq and phase lags —(f) of the system response at all temperatures. In Fig.ll, x 
and —(p are also plotted along with < W >. It is clear from the figure that < W >, Xq, and 
—<f) all peak almost at the same temperature. Interestingly, the peak value of —(f) is about 
0.827T. For this system —(f) satisfies peaking criterion of SR stated using the linear response 



theory 26) . However, magnitude wise the phase lag —<fi is off by about tt/2. 
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IV. DISCUSSION AND CONCLUSION 



A periodic potential system driven by a periodic applied field of small sub-threshold 
amplitude at a high frequency, close to the natural frequency of the periodic potential well 
bottoms, shows stochastic resonance. Here the average input energy per period of the 
field is considered as the quantifier of SR. The same quantifier had served SR correctly in 
bistable systems. Moreover, the probability distribution of the input energy exhibits similar 
qualitative behaviour at SR as it shows in the bistable systems. The linear response theory 
calculated frequency dependent mobility was found to show exactly similar behaviour as 
the input energy and was previously termed as a mere dynamical resonance. It was argued 
that the period of the drive was too small compared to the Kramers time of inter-well 
static potential barrier crossing, thus disqualifying the resonance behaviour from being a 
genuine stochastic resonance. However, we find that such an argument just does not hold 
because, in the dynamical situation, the inter-well transitions become quite numerous at 
the temperature where input energy peaks. Moreover, and most importantly, the intra-well 
trajectories show bistability and the obtained resonance is a result of transition between 
these dynamical states effectively helped by inter-well transitions. 

There has been two rivalling SR criteria involving the phase lag between the response 



and the applied field: one 



25] statin; 



phase lag equaling 7r/4 and the other 



that the phase lag shows inflection at SR with the 



26| that phase lag shows a peak at SR. We find that 
whereas the former criterion is only approximately satisfied in the uniform friction case, the 
latter is satisfied for the nonuniform friction case. Thus, these criteria seem to be true for 
specific systems and hence the phase lag cannot be taken as a universal quantifier such as 
the input energy for SR. 
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